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We examine electronic transport in a spin-blockaded double quantum dot. We show that by tuning 
the strength of the spin-orbit interaction the current flowing through the double dot exhibits a dip at 
zero magnetic field or a peak at a magnetic field for which the two-electron energy levels anticross. 
^^ ' This behaviour is due to the dependence of the singlet-triplet mixing on the field and spin-orbit 

amplitude. We derive approximate expressions for the current as a function of the amplitudes of 
the states involved in the transport. We also consider an alternative model that takes into account 
a finite number of nuclear spins and study the resulting coupled dynamics between electron and 



o 



JjH I nuclear spins. We show that if the spin ensemble is in a thermal state there are regular oscillations 

in the transient current followed by quasi-chaotic revivals akin to those seen in a thermal Jaynes- 
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Cummings model. 



INTRODUCTION 



' i _ A double quantum dot (DQD) can be used for the detailed investigation of spin interactions among electron spins 

Cj . and even between electron and nuclear spins. The interactions can be probed optically or electrically by monitoring 

I ■ the electrical current flowing through the DQD as a function of the energy offset between the two dots and the 

^ \ applied magnetic field. In a spin-blockaded DQD, the current through the DQD is large when the two electrons form 

a singlet state, whereas it is suppressed when the two electrons form a triplet statei The spin blockade mechanism 

is due to the Pauli principle and has been demonstrated in semiconductor heterostructure quantum dots^ as well 

as carbon nanotube dots.^ It has also been shown that the transient behaviour of the leakage current can provide 

valuable information about the interactions between electrons in the DQD and nuclear spins in the host material.^ In 

particular, the electrical transport process can lead to a coupled electron-nuclear dynamics, nuclear spin polarization, 

F^ ■ and hysteresis effects. ^^^^^ 

^ , In the spin blockade regime, the small leakage current increases when there is a process that leads to singlet-triplet 

O ■ hybridization (mixing). A non spin-conserving interdot tunnelling is one such process. This type of tunnelling may 

^ O I result from a spin-orbit interaction (SOI), which in some cases can be strong enough and thus has to be taken into 

account i^i^ The hyperfine interaction (HI) between dot electrons and nuclear spins can also lead to singlet-triplet 

hybridization. In a simplified approach, the nuclear spins create an effective magnetic field which acts on the electron 

^ . spins. This field can point in an arbitrary direction and mixes singlet and triplet states. 

On ■ The electrical current in a DQD system has been examined theoretically in the presence of strong SOI and in the 

^^ . regime where the coupling of the DQD to the leads corresponds to the largest rate in the system.^ In the first part of 

this work, we consider a DQD weakly coupled to the leads and examine the current in the resonant regime, i.e., when 

the lowest singlet and triplet energy levels are almost aligned. Unlike the approach followed in Ref. 7 we consider 

^J ■ explicitly the lowest one-electron states and derive rate equations that involve the transition rates between one- and 

two-electron states. The SOI is modelled with a non spin-conserving tunnel coupling amplitude tgo between the two 

dot^ii^ which couples the triplet states |T+), |T_) with singlet states. We consider weak SOI, that is tgo < ^c, where 

tc is the spin-conserving interdot tunnel coupling. The effect of a strong SOI on a spin-blockaded DQD has been 

investigated in Ref. with a more general SOI model in which all triplet states couple to singlet states. Further, we 

here give emphasis to the regime where the HI is weak enough or absent, which might be the case in carbon-based 

^ \ quantum dots.^^^^ The main aim of the first part of this work is to determine the current as a function of the amplitudes 

C3_' of the one- and two-electron states which participate in the transport cycle. Some approximate results which give 

valuable insight into the basic behavior of the current are derived. We show that, depending on the strength of the 

SOI, the current shows a dip at zero magnetic field or a peak when the lowest two-electron energy levels anticross. 

This behaviour occurs because when tgo is large the singlet-triplet mixing near zero field is much weaker compared to 

that at high field. This gives rise to a dip at zero field. However, when tgo is small the mixing is strong only near the 

anticrossing point leading to a peak in the current. 

In the second part of this work we focus mostly on the interplay between SOI and HI. To properly account for the 
HI, we employ a microscopic 'toy' model which takes into account a finite number of nuclear spins. To make this 
model tractable, only the two lowest singlet-triplet states are considered, and we treat the nuclear spins as a single 
large spin. We first consider the interplay between the SOI and HI, and its effect on the steady state transport and 
nuclear spin polarization. We find a sharp transition in the current and polarization as the SOI is increased, consistent 
with the topological phase transition investigated in Ref. [ll|. Second, we look at the transient dynamics induced by 



the HI alone, and find a strong oscillatory contribution depending on the hyperfine coupling strength and inversely 
proportional to the square root of the number of nuclear spins. 

II. ELECTRICAL CURRENT IN THE SPIN BLOCKADE REGIME 

A. Physical Model 

In this section the electrical current through the DQD is examined when the electron-nuclear spin dynamics are 
uncoupled. This section is concerned with the effect of the SOI on the electrical current while Hl-induced effects due 
to the coupled electron-nuclear dynamics are addressed in the next section. The DQD is modelled with the two-site 
Hamiltonian 

2 2 ^2 

Hbqd = i^c + Hso + i^hf + ^ SiTii + ^ ^ rii^riii + ^ ^ cr^- (1) 

Here H^ is the tunnel-coupling Hamiltonian that conserves spin and has the form 

He = -tc(4t^2t + cJ^C2;) + H.c, (2) 

and the Hamiltonian part due to the SOI that allows spin- flip has the form- 

^so = -tso{c[^C2l - c[^C2^) + H.C. (3) 

For the HI we assume the forn>i^ 

^hf = -^eMB(BN,lO-i + Bn,2^2). (4) 

The operator cj^ (qct) creates (destroys) an electron on dot i = 1,2 with spin a = {t,i} and orbital energy Si. The 
number operator is denoted by n^ = ^^ riia = cL^^t + 4i^*i- "^^^ tunnel coupling amplitude between the two dots is 
tc, the amplitude due to the SOI is tgo, the charging energy is [/, and A = g^ii^B is the Zeeman splitting due to the 
external magnetic field B in the z-direction. Here Bn,^ is the magnetic field in the zth dot due to the nuclear spins 
and (Ti are the Pauli operators. 

The quantum states which participate in the transport cycle in the spin blockade regime are the 2 lowest one- 
electron states and the 5 lowest two-electron states. For simplicity the three-electron states are neglected in this study 
because they do not change qualitatively the basic results. The one-electron eigenstates |j;n = 1) can be written in 
the general form 

|j;l)=a,|t,0) + /3,|0,t)+7.-U,0) + (^,|0,;), (5) 

with j = 1,2 and c]^cr|0) = \^->^)-> ^2^1^) = lO^cr). Here the eigenstates are ordered with increasing energy. In the 
spin-blockade regime U ^ t^ and further there is an energy offset between the two dots. In this work we choose for 
the on site energies 82 = Si — U/2 and define the energy detuning as ^ = ^(1, 1) — £^(0,2), where E{n^m) is the 
energy of the charge state which has n (m) electrons on dot i = 1 {i = 2). If H^o = and H^ = then ai = /3i = 
and 72 = ^2 = and the nonzero amplitudes satisfy Si ^ 71 and P2 ^0^2- When Hso 7^ the amplitudes Ofi, /3i, 72, 
S2 are in general nonzero and satisfy /3i <C o^i, S2 <C 72. 

Neglecting double occupation on dot 1, a two-electron eigenstate \j;n = 2) with j = 3, ...7 has the general form 

|j;2) = a,|t,t) + ^,|t4) + c,|;,t)+c^,|;4) + e,|0,n). (6) 

Here |T_) = | ;,;) = c\^cl^\0), |T+) = | t,t) = c\^cl^\0), \a,a^) = clcljO), and \So2) = |0,n) = 4t4j0)- The 
amplitudes of the various components depend on the strengths of the SOI and HI as well as the Zeeman splitting 
and detuning. When H^o = and H^ = the two-electron eigenstates correspond to the triplet states |T+), 
|T_), |To) = (I ti) + I it))/A/2 and the two singlet states, which consist of the components |6'o2) = |0,ti) and 
1^11) = (I ti) ~ I it))/>/2- The effect of the J^so is to couple the |T+), |T_) states to singlet components. As shown 
below the coupling strength increases with tgo and, for a fixed detuning, is sensitive to the Zeeman splitting. 

To calculate the electrical current flowing through the DQD we employ a density matrix approach within the Born 
and Markov approximations.^^ The internal parameters of the DQD and the chemical potentials of the two leads are 
adjusted to the spin-blockade regime. This regime can be identified from the fact that when H^o = and JT^hf = the 
three triplet states |T+), |T_), and |To) are equally and almost fully populated (~ 1/3), provided spin decoherence is 
ignored, and the current as a function of the source-drain bias is suppressed. 
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FIG. 1: Electrical current, /, as a function of the Zeeman splitting, A, and energy detuning, 6, when there is no hyperfine 
interaction. The SOI amplitude is tso = O.Oltc, where tc is the spin-conserving interdot tunnel coupling. 



Results 



When tso is nonzero the |T+), |T_) states couple to singlet states and this coupling has a direct effect on the current. 
To demonstrate this effect we show in Fig. [1] the current as a function of the Zeeman splitting A and energy detuning 
S between the two dots when there is no HI. The coupling of the states is strong near the anti-crossing points leading 
to an increase in the current (see also below). As a result the curves of high current map-out the points where the 
energy levels of the quasi singlet and triplet states anticross. When t^o = the leakage current is approximately 
constant and no high current curves occur. 

To understand the role of the SOI we show in Fig. [2] the energy diagram of the two-electron states as well as the 
current as a function of the Zeeman splitting A, for a fixed energy detuning S and different SOI amplitudes with 
Hm = 0. We concentrate on the regime tso < ^c and choose (5 < 0, which is experimentally the most interesting 
case for spin-qubit application^ i^^i^^ , because the spin pair can be described by an effective Heisenberg model. The 
resonant current that is defined as the current at the anticrossing point increases with tso and the same occurs for the 
asymptotic current that is defined as the current at a high magnetic field. Therefore, when tso is large enough the 
asymptotic current becomes approximately equal to the resonant current and thus the peak cannot be distinguished. 
The same pattern occurs when the direction of the magnetic field is reversed, thus the current as a function of the 
Zeeman splitting shows a dip at A = 0. To a good approximation this pattern is independent of the detuning, 
provided tso is large, and consequently the anticrossing points of the energy diagram cannot be probed. 

To quantify the above results we analyze the rate equations and calculate analytically the transition rates. Then 
we derive the steady-state current for some interesting limits, such as for example the current at the singlet-triplet 
anticrossing point, as well as foiB^O and B high. We are interested in determining the current in the steady-state 
for a DQD that is weakly coupled to the leads. In this regime we can consider only the diagonal elements of the 
density matrix^^ and the dynamics of the system is described by the rate equations 






dt 



/ ^ Pra-t^mn-) 



(7) 



where the diagonal elements of the density matrix are denoted by pn = pnn and the normalization condition is 
Sn=i pn = ^' The transition rate from an eigenstate \n) of the DQD with eigenenergy En to an eigenstate \m) with 
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FIG. 2: (Color online) The top panel shows the energies, E^ of the four lowest two-electron states as a function of the Zeeman 
splitting, A, for an energy detuning 8 — — 25tc and SOI amplitude tso = O.ltc, where tc is the spin-conserving interdot tunnel 
coupling. The notation of the states is given in the main text in Eq. JB]). The bottom panel shows the electrical current, /, as 
a function of the Zeeman splitting A for an energy detuning 8 — — 25tc. 



eigenenergy Em is 



Rum = r^[|(n|ci,|m)p/L(^mn) + |(m|ci,|n)p/Z(^nm) 
+ \{n\c2a\m)\''fR{E^n) + \{m\c2a\n)\^ fR{Enn,% 



(8) 

where V = 27r\tr\'^Dr/h. The tunnelling amplitude between dot and lead is t^, fi{Emn) is the Fermi-Dirac distribution 
function at chemical potential /i^, with Emn = Em — En and /~ = 1 — /. Also, Dr is the density of states for the 
leads, which we assume to be constant and equal for both leads. 

The operator for the electrical current, for example, for the right lead is 



i = T^^tr{c\^dka - dl^C2cj), 



(9) 



k(T 



where dka denotes a lead operator. Tracing out the leads we derive the following expression for the average current 

/ = eF ^ Pn[\{n\c2^m)\^fR{Emn) - \{m\c2M? fA^nm)]- (10) 

nmcT 

Starting with the rate equations and calculating the transition rates it can be readily derived that the absolute value 
of the current for H\^i = is / = e.^^^^{Rni + Rn2)pn- The simplest regime is when tgo = and 5 = 0. Here, the 
occupations of the triplet states satisfy Pa = p5 = pe = Pt- From the steady-state condition, dpn/dt = 0, it can be 
derived that 
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FIG. 3: (Color online) The figure shows the important quantities which determine the transition rates through the parameters 
Fo, ^0 (top panel) and Fa, Ea (bottom panel) defined in Eqs. (|14p , (|15p and Eqs. (|2Qp , (|2ip respectively. The energy detuning 
is ^ = — 25tc and the Zeeman splitting is chosen in the low field regime A = 0.02tc (top) and asymptotic regime A = 0.7tc 
(bottom). 



For weakly coupled dots the second term is typically negligible and pt ^ 1/3. It is easy to prove that this approxi- 
mation is excellent when tc is small and ^ is large. In this regime the leakage current is 
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For tgo 7^ we calculate analytically the transition rates and derive that for low magnetic fields the current is given 
approximately by the expression 
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FIG. 4: (Color online) The figure shows the absolute amplitudes 63, 64 and ee which express the mixing of the triplets |T+) and 
|T_) with the |*S'o2) state. The Zeeman splitting is chosen in the low field regime A = 0.02tc and asymptotic regime A = 0.7tc. 



The transition rates i?i5, R25 involve the one-electron state |1; 1), |2; 1) respectively, and the triplet state |To) = |5; 2). 
The terms which are proportional to R12, R25 do not affect the physics we examine here, so for this reason they 
are not given explicitly. Also, as seen in Fig. [21 in the small tso regime that we are interested in /q is to a good 
approximation independent of tso- An approximate expression for the resonant current is 



with the parameters 
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(18) 



The resonant current corresponds to the magnetic field B^ for which the lowest quasi singlet and triplet states anticross 
and it is well-defined for t^o <C tc- The asymptotic current that corresponds to a high B is given approximately by 
the expression 



with 
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The asymptotic current is defined at a high magnetic field where the current varies slowly with B. 

Equation (p!3|) can be used to estimate the width of the current peak that occurs at A^ = g^jj^BBr- Specifically, if 
we denote by A^ — Ai the half width at half maximum of the peak, then Ai satisfies the relation /o(^i) = ^r/2. In 
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FIG. 5: (Color online) Electrical current, /, as a function of the Zeeman splitting, A, for the energy detuning S 
different combinations of SOI and HI strengths. Here A^ = ^e/iscTN. 



-25tc and 



the same way the width of the A = minimum can be estimated in the regime la ^ Ir- If Ai corresponds to half of 
the width then to a good approximation it has to satisfy the relation /o(Ai) = /a/2. 

The variation of the current as a function of the magnetic field and SOI strength is mainly due to the change of 
the first terms in Ei and Fi. For example, if tgo <C tc, then as can be seen in Fig. [3] only the quantities (71(^4)^, 
(a2<^6)^7 cind (71(^3)^ which contribute to Ei and Fi are important. Because these quantities are approximately equal 
we derive that Ea ^ Eq and Fa ^ Fq^ thus la ^ lo- In the same way, at Br the corresponding amplitudes in Eqs. ([TTj) 
and ([iHj) lead to £^0 > Er and Fq > Fr. Thus, Ir > Iq and a peak is formed at Br. On the other hand, in the range 
4o < tc the quantities ((^163)^, (726^3 + (^2^3)'^^ and (/32e6)^, (o^iae — ^lee)^ increase drastically (Fig. [3]). As a result, 
Ea and Fa decrease significantly, whereas £^0 and Fq do not change a lot. Therefore, for intermediate or large tgo the 
asymptotic current la is much larger than the current at very low fields and it approaches the current Ir. Eventually 
as tso increases the current shows a dip for A = 0. 

Inspection of the various amplitudes involved in the transition rates demonstrates that the important amplitudes 
in order to understand the current are the 64 (or 63 when B is high) and cq. These amplitudes express the mixing 
of the |T+) and |T_) states with the |5'o2) component. These amplitudes are responsible for the behaviour of the 
current when tso is included in the model. As shown in Fig. IH when tso is not too small the mixing of the |T+), |T_) 
states with the 1502) component is stronger in the asymptotic regime than for B ^ 0. For this reason the current 
is more sensitive to ^o when B is high but shows only a small variation with tso for B ^ 0. Furthermore, as B 
increases (in the asymptotic regime) then for a fixed tso and negative detuning the amplitude cq for the |T+) state 
increases, whereas the amplitude 63 for the |T_) state decreases. This behaviour can be understood by noticing that 
the corresponding anticrossing points move along the detuning axis in opposite directions. In Eq. ([19]), Ea + Fa is to 
a good approximation constant and the current la remains constant as A increases. 

To take into account the hyperfine interaction, we follow a standard approach and treat the nuclear magnetic fields 
BN,i in the two dots as quasi-static classical variables which take random values^ In this case the electron and 
nuclear spin dynamics are uncoupled. In the next section we develop a model to look at the coupled dynamics. The 
distribution for each random static field is Gaussian with spread ctn- The electrical current is computed as the average 
over different random field configurations^ This is a good approximation when the nuclear dynamics is much slower 
than the electron dynamics. In Fig.[5]we plot the current versus Zeeman splitting. When tso = 0, a peak is formed at 
the singlet-triplet anticrossing point due to the mixing caused now by the HI. When both tso and (Jn are nonzero the 
resonant and the asymptotic current increase. However, for A = the current is determined by the HI, consistent 
with the results in Ref. 7. As the spread cr-^ increases, the HI results in a peak at A = 0, though the A = current 
may have a more complicated form when (Jn is large. In the small <jn regime that we are interested in and for tgo 7^ 0, 
the A = current is typically smaller than the asymptotic current. Our numerical calculations confirm that these 
trends occur for different choices of tso provided that gefJ^BCfN "C tso- 



III. NUCLEAR SPIN POLARIZATION AND TRANSIENT DYNAMICS 

A. Physical model 

The model employed in the previous section provides insight into the current, but it does not capture the coupled 
electron-nuclear spin dynamics. Thus, for example, the nuclear spin polarization as a result of the transported 
electrons through the DQD cannot be examined. In this section we look at the DQD system from a different 
perspective. Specifically, we use an idealized model to study the coupled electron-nuclear spin dynamics and how 
this affects the transient behaviour of the current. We will use this model to show two things. First, that the 
presence of even a weak spin-orbit coupling can prevent nuclear spin polarization. Second, if the nuclear spin state is 
completely thermalized several interesting features arise in the transient electron current; regular beating followed by 
quasi-chaotic oscillations. 

A true model of the 2^ states which describe the N nuclear spins is computationally intractable, but some insight 
can be gained from, e.g., treating the nuclear spins as a "giant" spin.^^^^^^^ In addition to reduce the complexity 
even further, we restrict ourselves to a small subspace of the two-electron Hilbert space, spanned by the states |T+), 
IS'ii), and |5'o2)- This is a reasonable approximation under an appropriate bias, i.e., when only the state |T+) is in 
the voltage bias window (and neglecting tunneling, from the reservoirs, into superpositions of the singlet states). In 
addition we assume that the rate of tunneling from the left lead to the dot is large, and that we are at the anti- 
crossing point of the singlet-triplet states.^ This reduction of the state-space does also neglect the occupation of the 
|0,^) state, and trapping in other single-electron-occupation states. As a test, we included the single-electron state, 
|0, 1), in an extended model, but it had little impact on the results we present here. Thus hereafter it will be neglected. 
In addition, another interesting regime to investigate would be to assume a larger bias and include all three triplet 
states. However, since all the interesting coupled electron-nuclear effects we discuss here are mediated by the |T+) 
and l^*!!) dynamics, such a regime may only reduce the visibility of these effects. 

We model the interaction of these three states with a non-equilibrium master equation. This model allows a fiow 
of electrons through the DQD, and thus we can estimate properties like the total polarization of the nuclear spin, 
and the transient dynamics of the coupled electron-nuclear spin system. This is in contrast to treating the spins as a 
frozen spin bath,^^ as a semi-classical degree of freedom, ^^^— or as a non-Markovian environment ^21 

The total Hamiltonian of the system is then written as 

F = iJe + ffso + Hm - A|r+)(T+| - S\So2){So2\, (22) 

where He describes the couphng between the singlet states 

H, = mSn){So2\ + \So2){Sii\], (23) 

and Hso is the Hamiltonian part due to the SOI, which couples singlet-triplet states, and has the form 

Hso = tso[\T+){So2\ + |^02)(T+|]. (24) 



The hyperfine interaction is 
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T+){T+\I^ + ^{a'^yi'l+a^^hl' 



(25) 



Here S is the energy detuning and A = ge/j^BB is the Zeeman splitting caused by the external magnetic field B in the 
z-direction, tgo is the spin-orbit coupling, and tc is the coherent tunnelling amplitude for the singlets. In JT^hf we use 
an effective-spin notation so that a\_^ = |T+)(S'ii|. The coupling g is the nuclear hyperfine coupling term, g = A/N^ 
where N is the number of nuclear spins and A is in the range of 80 /ieV. 

For simplicity, we assume the nuclear spins are spin- 1/2, and the hyperfine coupling strength is uniform and 
homogeneous. This implies equal-size dots, and an equal nuclear hyperfine coupling on each site. We are working 
in a rotating frame and the sign difference between couplings on the left and right dots, due to the antisymmetry of 
the wavefunction, is omitted. ^^ Thus the nuclear states are in fact the difference between nuclear spin states in the 
left dot and the right dot. This allows us to use the giant-spin approximation Ji = J]]^ /^^, i = z, +, — . For a large 
thermal state one should really describe the spi n sys tem as a distribution over giant spins of differing sizes. ^^ Here 
we consider a single giant spin of length J = ^/N/2, which may be valid if the distribution of spin sizes is strongly 
peaked (see Ref. Tt for a rigorous discussion of this assumption). In the final section we discuss the possible effects of 
a broader distribution. 



To account for electron transport we include a Lindblad term, 

r 

Li[p\ = 2 [^^+^^- ~ /^-/^+P - PM-M+] ^ (26) 

where the electron transport operator is /i_ = |5'o2)(^+|, and we have omitted the vacuum state which is a good 
approximation when the rate of tunnelling- in from the left lead is large. Thus, in the results that we show in the 
following figures, we solve the master equation for the combined electron- nuclear spin density matrix p, 

f = 4[^,p]+i:.b]. (27) 

Because of the small Zeeman splitting of the nuclear spins, we impose an effective infinite temperature for the nuclear 
system initial state. Nuclear spin dephasing and thermalization are not explicitly included here as they can, in 
principle, break the large spin approximation, and cannot be written in terms of the Jz operators alone. Fortunately, 
the nuclear spins are typically weakly interacting with each other and the environment. We assume that the current 
is dominated by transport into the right reservoir, and is defined as 

I{t)=eTT[Jnp{t)], (28) 

where the superoperator is the jump operator Jn[-] = r/i+ • /i_. For the polarization of the nuclear spins we calculate 
the expectation value of Jz for the large spin. 

B. Nuclear spin polarization 

Hysteresis in the current measurement, as one sweeps the external magnetic field through the singlet-triplet level 
crossing, is a sign of nuclear spin polarization. In vertical dots, large polarizations (> 40%) have been achieved^"— 
In lateral dots, the polarizations are significantly smaller, perhaps due to either the asymmetry in the dots, and thus 
in the nuclear hyperfine interaction, or dark states. ^^ However, even for small levels of nuclear spin polarization, large 
hysteresis has been observed. Further, in several recent studies it has been shown theoretically that if the spin-orbit 
coupling is above a certain threshold, relative to the nuclear hyperfine coupling, no polarization of the nuclear spin 
occurs. ^^ If it is below that threshold, the nuclear system becomes polarized due to the spin- flip process during the 
electron transport. At some critical value between these two regimes, long-lived dark states can occur, alongside a 
topological phase transition. 

To investigate this phenomenon in our model we plot, in Fig.[6l both the steady-state current / and the normalized 
nuclear polarization Jz/J as a function of tso for a given nuclear hyperfine coupling A = 0.1 meV. We use dot 
parameters which put us at the singlet-triplet resonance point, and also omit the first term (the Overhauser term) in 
Eq. (|25]). We see that for 



Uo < 25A/4V7V, (29) 

the nuclear spin is strongly polarized by the electron transport process, and the current flow is low. Conversely, as 

tso > 2bA/4^/N, (30) 

the spin-orbit mediated transport route becomes dominant, and the nuclear spin is no longer maximally polarized. 
The transition between these two regimes is consistent with the sharp change observed in the "topological phase 
transition" ^^. The factor of 25 arises from the amplitudes Co, 2 and Ci,i of the combined effective singlet state used in 
that treatment \S) = Co,2|'S'(0, 2)) + d,i|6'(l, 1)). Inspection of the bare eigenstates of the electron spin Hamiltonian 
shows that this factor is, in general, 

|(Cl,l/Co,2)| = |(<^ - \/4tfT#)/2ie|. (31) 

Thus the quantum dot system parameters can in practice also be tuned to sweep the axis of Fig. [6l In the bottom 
half of Fig. [6] we show the absolute value of the matrix elements (in the basis of the m G { — J^ J} eigenstates of Jz) 
of the steady-state nuclear spin density matrix for three choices of spin-orbit coupling: zero, at the "critical point" 
and far above the critical point. One can easily see that at the critical point there are significant non-zero coherences 
in the nuclear spin state, which may be related to the presence of long-lived "dark states". ^^ Our results indicate 
that precursors to this "topological phase transition" exist even in the presence of a full transport cycle and noisy 
environmental effects, akin to the persistence of quantum phase transitions in non-equilibrium systems J^ 
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FIG. 6: (Color online) Top: Steady-state current, /, (solid lines) and nuclear spin polarization, J^, (dashed lines) as a function 
of tso. Three different sizes of large spin are shown, J = 10 (blue), J — 20 (red), and J = 40 (black), corresponding to A^ = 200, 
400 and 3800 nuclear spins respectively. The transition occurs for tso > 25A/4A/iV, where A is the maximum value of the 
nuclear hyperfine coupling g = A/N. Because of the inverse scaling of the transition point, as a function of N, these results 
suggest that a relatively small spin-orbit coupling can suppress the polarization of the nuclear spin. Bottom: matrix elements 
(in the m basis of the Jz operator) of the steady-state reduced density matrix of the nuclear spin, for J = 20. At the critical 
point the steady- state contains significant non-zero coherences. 



C. Transient dynamics 

Some evidence^^ indicates that the low dephasing and extremely long relaxation time of the nuclear spins, 
combined with the fast stochastic electron transport process inducing nuclear polarization, leads to a long time 
instability phenomenon and fluctuations in the nuclear spin state. ^^^^^ This is typically observed in the long-time 
beating seen in the current. ^^^^^ This is almost certainly a semi-classical eSect^^^^^^ though some evidence 
suggests that coherence within the nuclear spins can survive for millisecond time scales^ To gain some insight on 
what might be observed on shorter time scales, we now examine the transient behaviour of the electron transport 
induced by the nuclear spin. In the following we entirely neglect the spin-orbit coupling. 

We assume that at some initial time t = the electronic system is prepared in the |T+) state and the nuclear system 
is in an initially maximally- mixed state 



P(0) 



2J+1^-^ 



\m){m\®\T+){T+ 



(32) 



We then solve the dynamics given by this initial state in the equation of motion. As we will show below, we find that 
this produces dynamics akin to the Jaynes-Cummings Hamiltonian from quantum optics with an initial thermal (or 
chaotic) cavity state. ^^'^^ This is a well-studied system, with so-called "quasi-chaotic revivals". In the optical case, 

27rh/Xy/E, 



studies have shown that one observes an initial sharp peak'^^ in the atomic state on a time-scale t 

where A is the field-atom coupling in that model, and n is the initial thermal occupation of the field given by the 



^(bosonic) 
P 
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FIG. 7: Transient current, /(t), with an initially maximally- mixed nuclear spin state for (a) J — 10, and (b) J — 40. The large 
number of commensurate frequencies around m — cause the regular oscillations with frequency A/4v]V, before the onset of 
quasi-chaotic dynamics. 



Bose-Einstein distribution. This is typically followed by the onset of the so-called "quasi-chaotic" oscillations at a 
time that scales with the thermal occupation. 

In Fig. [7] we show the transient electron current for two sizes of the nuclear spin system. Again we omit the 
Overhauser term. In the figure, one can see the clear onset of quasi-chaotic beating and a scaling of the onset of this 
beating with the size of the large spin J. Since we are now dealing with an infinite-temperature large spin, the role 
played by n in the optical case is now played by J in the large-spin case. Unlike in the optical case, however, the 
collapse of an initial peak is not seen here. Instead, there is a transition from large steady oscillations to quasi-chaotic 
ones. 

To understand the presence of the large steady oscillations we can make two observations. Firstly, for low thermal 
occupation n <^ J ^ and large J, the Holstein-Primakoff transformation tells us that the dynamics would coincide with 
that of the bosonic Jaynes-Cummings model, but with renormalized coupling A ~ Aj^^plN . As n oc J^ the analogy 
should break-down. Secondly, the large regular oscillations can be understood by explicit diagonalization of the iJhf 
Hamiltonian (again omitting the first, Overhauser, term), which gives the following expression for the occupation 



probability of the |T+) state, assuming that the initial state is |'0(O)) = |T+) ^J^ 



'Cm\m), 



\{T^mm' 



J 



\c„ 



|2 2 

COS 



E+t 



(33) 



where 



E. 



(m) ^ Av/J(J + l)-m(m + l) 

2V2N 



(34) 



In the bosonic case, this is an infinite sum whose frequencies scale as vrH-T. One can understand the large regular 
oscillations that occur in the large spin case by considering the contributions to the sum around m = 0. This 
region contributes a large number of terms to the sum, with similar frequencies E^^ ' jh ^ Aj^h^fN ^ which are 
commensurate at small times. This gives rise to the frequency of the early-time oscillations in Fig. [71 

The possibility to observe both the early-time oscillations and the quasi-chaotic oscillations in experiment is in- 
triguing. One can consider that as N is increased the period of the oscillations increases. For N ex 10^-10^ this can 
reach the order of 100 micro-seconds. Ultimately, the observation of these oscillations is limited by two effects. First, 
electron spin dephasing will affect these dynamics. We found that, by including such effects in our model, if these 
dephasing rates are much larger than the frequency E^ the oscillations become damped. However, the oscillations 
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are not strongly affected by dephasing or decoherence in the charge degr ee of freedom. Secondly, in reality, the width 
of the distribution of large spins in the thermal ensemble around \/N/2 will also cause dephasing, and is the primary 
ca use o f the electron spin dephasing to begin with. One can consider that large spins in this distribution close to 
\/N/2 will contribute "in phase" to the early-time oscillations and those far away will induce additional dephasing. 
Thus, inevitably the oscillations shown in Fig. [71 will decay if this distribution is broad. Finally, dephasing and 
rethermalization of the nuclear spin states will also affect the dynamics if these rates become comparable to Ej^. 

In addition, the behavior we show in the dynamics of the current occurs on a relatively short time scale and strongly 
depends on the initial state. We also emphasize that the meaning of this short-time dynamical current we plot in the 
figures is the following: it is the ensemble average of detecting the current of a single electron entering the reservoir 
based on a specific initial condition. In other words, in a real experiment, one would have to repeatedly prepare the 
nuclear-spin and dot system in the same initial state, and measure the resulting single electron transport events to 
eventually collate the data shown in the figure. Obviously this is experimentally challenging. A more feasible and 
natural approach would be to detect transient behavior in the high-frequency current-noise spectrum. However, given 
that such transients are "around the steady state" (a state which may include significant nuclear spin polarization), 
some of the features which rely on the nuclear spin state being in a maximally-mixed state may become less visible. 

In practice, it maybe more feasible to consider a closed system, i.e., disconnected from electronic reservoirs, and 
measured by, e.g., a charge detector. In this case there will be no dynamical nuclear polarization, which is advantageous 
for observing the features which rely on the maximally mixed state. Also, at this stage, it is not clear if there is any 
connection between the oscillations we observe here and those seen in experiments. Even for large N the time scales 
still differ greatly. Finally, an alternative system to investigate this phenomenon could be quantum dots in carbon 
nanotubes, or with superconducting qubits/wave guides^^ coupled to ensembles of spins in diamond. 

IV. CONCLUSIONS 

This work investigated electronic transport in a double quantum dot for weak spin-orbit interaction, namely, when 
the spin-orbit amplitude is smaller than the interdot coupling. The electrical current was calculated numerically from 
rate equations. We derived simple approximate expressions for the current as a function of the amplitudes of the one- 
and two-electron states which participate in the transport cycle. We found that when the SOI is small the current 
shows a peak at a magnetic field for which the lowest two-electron energy levels anticross, whereas when the SOI is 
large a dip is formed at zero magnetic field. Numerical calculations showed that in a double dot system with a small 
hyperfine interaction these characteristics remain valid. 

We also considered a model which includes dynamical behavior of the nuclear spin, and investigated the coupled 
dynamics between electron and nuclear spins. We found that the conflict between the spin-orbit and nuclear hyperfine 
couplings results in a transition point between no polarization and large polarization of the nuclear spin ensemble. 
We also considered the transient dynamics, where the nuclear spin ensemble is initially prepared in a highly thermal 
state. The unique characteristics of the large effective-spin model used to describe the nuclear spin ensemble induces 
dynamics in the current which depends on the fundamental nuclear hyperfine coupling and the ensemble size. 
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